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ABSTRACT 



This paper deals with the refraction of long 
gravity waves about a circular island with uniform 
slope below mean water level. The Fermat principle is 
used as a basis to derive an analytical solution for 
rhe ray paths. A refraction analysis is performed to 
determine the refraction factor, phase lag, and wave 
height at the island shore. The relative power distri- 
bution of rays in the far field is determined by con- 
sidering emergent reflected and refracted waves from the 
island . 
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CHAPTER I 
INTRODUCTION 

lo Statement of the Problem 

The purpose of this paper is to study the influence 
on long gravity waves of an island with circular bottom 
contours and radially constant slope- Using the prin- 
ciples of wave ray theory, a refraction analysis is 
carried out to determine the modification of initially 
plane waves as they move into shallow water about the 
island. The refraction and direct shoaling effects 
upon incoming long water waves are considered. The re- 
flected and refracted rays emerging from the island 
are considered in order to study the divergent wave 
energy in the far field of the island. This refraction 
study is a contribution in part to a more exacting treat- 
ment of the subject by Vastano (1S66), involving a 
numerical diffraction technique. 
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2. Historical review 



For centuries man has concerned himself with the 
effects of islands on surface waves. The early civiliza- 
tions of Micronesians and Polynesians utilized wave pat- 
terns to mark fishing grounds and aid in navigation 
among the islands. Winkler (1901) recorded some of their 
techniques through interviews with native chiefs * Kriege 
(1943) further illustrates this fact and states: 

Every chief and sea pilot possessed 
elaborate charts based upon his own exper- 
ience and on knowledge handed down or gain- 
ed from others. These sea charts are made 
with thin strips of the midrib of the leaf- 
lets of the coconut arranged on a frame 
usually rectangular in shape. The knowledge 
they record is indicated by the arrange- 
ment of the leaf strips relative to one 
another and by the forms given to them 
by bending and crossing. Curved strips in- 
dicate the altered direction taken by ocean 
swells when deflected by the presence of 
an island; their interactions are nodes 
where these meet and tend to produce a confused 
sea. 

More recent interest on the specific subject of wave 
refraction was stimulated by the needs brought about by 
'World War II amphibious operations. As a wave approaches 
shallow water, its velocity begins to decrease appre- 
ciably at a depth of about one-half the wave length and 
the wave height is influenced by the convergence and div- 
ergence of energy along the shore. As the wave bends 
around the island, "wave shadow" areas may be formed in 



but often times not to the extent 



the island "wave lee," 
that might be expected. This fact was painfully realized 
during the invasion of Sicily when the methods of Sverdrup 
and Kunk (1947) were used to forecast 'wave conditions 
prior to the landing. Despite the use of these acceptable 
methods, the wave heights in the lee of the island were 
greatly underestimated thus severely hampering the landing 
operation „ 

The graphical means of constructing wave refraction 
diagrams was first used by O'Brien (1942). Application 
to the complex ocean bottoms of nature was made by 
Johnson, O’Brien, and Isaacs (1948) and Arthur, Munk, and 
Isaacs (1952). The wave refraction pattern for any bottom 
topography assumes such that the waves change direction 
and are bent in such a manner that they tend to assume 
the slope of the depth contours. Aerial photographs 
(Munk and Arthur, 1952) qualitatively bear this out. 

It was learned that the effect of bottom features upon 
v/aves depends not upon the absolute depth of water, but 
upon the depth relative to the wave length. Graphical 
construction methods by hand are both time consuming and 
tedious. Griswold (1963) has computer programmed the 
graphical construction of wave rays and Momoi (1964) has 
expanded on Griswold's program to include the effects of 
v/ave refraction. 



In some special cases, analytical solutions to the 
wave refraction problem have been developed to improve 
accuracy which can be critical in the investigation of the 
energy distribution especially on the wave lee side of is- 
lands o Analytical and experimental studies of refraction 
in the region of a caustic have been made by Pierson (1951) 
A method of conformal transformation was applied by 
Pocinski (1950) « The actual cases treated are restricted 
by the mathematical complications that are introduced by 
a particular conformal transformation and by how well 
the transform simulates a beach contour. Arthur (1950) 
used Fermat’s principle as a basis for an analogy between 
the refraction and the minimal flight path problem. The 
effects of current and depth distribution are considered. 
Johnson (1947) also discusses the refraction of ocean 
surface waves by current for deep water v/aves, 

Nunk and Arthur (1952) applied ray theory to calcu- 
late the wave intensity along a refracted ray. This pro- 
cedure has certain advantages over the usual technique 
of computing intensity from measurements of the distance 
between rays; namely in the areas of extreme convergence 
or divergence. 

Analytical solutions for the refraction of waves 
around an island with concentric ci cular bottom contours 
have been obtained by Arthur (1945) through the application 
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of Fermat 1 
important 
the "wave 
refraction 
direction. 



s principle, Arthur (1951) determined n’.ac the 
effects in the penetration of wave energy into 
lee" of an island are generally the resu.t or 
by underwater topography and variability in 
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3 . General Description of the Islanc 

The island bottom configure cion selected for the 
present study is axially symmetrical (i.e., the depth 
contours are concentric circles.) At 5 km. radius from 
the island axis the topography drops vertically from 
mean sea level to a depth of 0.5 km., then assumes a slope 
of 1:10 out to 41 km. radius, beyond which the depth 
remains 4.1 km. (see Figure 1). The periods of the 
waves selected are 2, 4, 8, and 12 minutes to provide a 
variation in wave lengths within the tsunami class. 

There was little attempt here to simulate the 
detailed underwater topographical features of any part- 
icular island. However, a computer program for the 
study of the interaction of monochromatic waves with an 
island of irregular shape is being developed by Reid 
and Vastano (1966). This program permits depth variation 
and the use of a special orthogonal system. 

It is interesting to note that the selected island 
configuration lends itself to a fairly simple mathemati- 
cal treatment. However, if a case were chosen where the 
extended slope line intersects the mean sea level at 
other than the island vertex, then a somewhat more for- 
midable problem presents itself with regard to wave 
amplitude determination and is worthy of further in- 
vestigation „ 
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4. Assumptions 

In water of constant depth the wave celerity, C, 
for sinusoidal gravity waves o_ small amplitude, neglect- 
ing surface tension, is given by the classical relation 
(Lamb, 1932) 

(1) = (g/lc) tan'n kh , 

where h is the mean depth of water, k = 2 77/L the wave 
number, and g the acceleration due to gravity . A wave 
schematic is shown in Figure 2, where the wave height h is 
equal to twice the amplitude, ^-. e . the vertical distance 
between two successive crests, and the period T is the 
time interval between the passage of two consecutive 
crests. For any wave greater than 0.1 feet in length, 
the effect of surface tension on wave speed can be ignored 
(Wiegel, 1964). The wave speed is related to wave length 
and period by 

(2) L = CT o 

There are two limiting cases of (1) which are 
appropriate to the so called "deep" and "shallow" water 
waves respectively. For deep water , h/L is very large 
and tanh kh tends toward unity. Under this condition (1) 
and (2) yield 



a 



2 

2 S- 9 1 • 

c = wr‘ \wrj ■ 

On the other hand, under the condition h/L — 1/20, (1) 

becomes to a good approximation 

C 2 = gh , 

There is, of course, no sharp distinction becween 
"deep" and "shallow" water,. The effect of depth of water 
on wave characteristics is gradual and waves in any 
finite depth of water are affected by that depth.. It is 
interesting to note that waves classified in the shorter 
period category are not refracted until near shore and 
shallov; water, while long period waves are refracted 
further offshore in deeper water. Tsunamis, for example, 
of periods 2 minutes to 1 hour would constantly undergo re- 
fraction effects in the open ocean. In a tsunami study, 
only waves that travel at a speed of (gh)^ are considered. 
As waves propagate into regions of changing depth, 
we consider that T remains unchanged while H and L vary. 
Given an initially plane wave of a given period and 
wave height moving towards shallow water, the wave height 
in general can vary according to the following factors; 

(a) Direct shoaling effects 

(b) Wave refraction effects 

(c) Reflection 
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(d) Diffraction 

(e) Energy lost oy fricaion and percolacicr. 

In this study, only the first two factors are considered 
over the sloping bathymetry, however, reflection is taken 
into account at the shore. The energy dispersion due to 
diffraction occurs in areas of appreciable convergence 

or divergence of energy and ir. shallow water zones such 
as behind breakwaters . The diffraction tends to spread 
the wave energy into potential shadow zones. Friction 
losses on slopes greater than 0.01 are considered neg- 
ligible (Walsh, Reid, and Bader, 1962). Percolation is 
related to the bottom characteristics and is considered 
negligible with respect to friction (Reid and 
Bretschneider , 1953). In addition, external wind influences 
are neglected. Xajiura (1964) discusses the partial 
reflection of water waves passing over a bottom of 
variable depth. Green’s formula provides a very good 
estimate of the transmitted wave amplitude even if 
partial reflection exists on a sloping bottom (Reid, 

1957) . 

In the absence of energy losses by friction and 
percolation or dispersion by. diffraction and reflection, 
the wave energy transmitted between orthogonal s (wave 
rays) remains constant. 
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CHAPTER II 
THEORETICAL ASPECTS OF THE PrtCuLLh 
1. Wave Refraction 

Gravity wave refraction is analagous to the bending 
of light rays in geometrical optics. Since the wave 
speed decreases with decreasing depth (C = (gh) 2 ) , rhe 
wave rays are refracted cowards regions of lesser water 
depth. 

This can be put in quantitative terms by the relation 



(1.1) 


&P _ 1 DC. 

ds C 0n 


where 


is the angle of the wave ray measured counter- 



clockwise from a fixed reference (the x-axis), s is arc 
distance on the wave ray and n is distance measured 
parallel to the wave crest and to the left of the ray. 
The above relation can be shown to be consistent with 



Fermat ' s 


principle (hunk and Arthur, 1952). In che spec- 


ial case 


of straight and parallel depth contours, one 


can take 


the x-axis normal to the contours such that 


C = C(x) 


alone. In this case (1.1) reduces to 

cos j3 i sin s 42 

' dx C r ox 


which is 


readily integrated and yields Snell's law- 
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Cl. 2) 



sin 



3 



C 



cons cent 



along the wave ray. Note however , that this is a _pecia 
case of the more general relation (1.1). 

In the general case, graphical methods are often 
referred to in describing the refraction pattern. The 
two methods most frequently employed are: first, the 
crest method, where successive positions of the wave 
crest are drawn by plotting the wave advance from point 
to point along the crest; second, the ray method, where 
each orthogonal is plotted directly by determining its 
shoreward deflection as it crosses successive bottom 
contours . Johnson, O'Brien, and Isaacs (194S) provide 
explicit instructions for this approach. In both cases 
a detailed large scale chart of the bottom topography 
is essential. For the advantages ar.d disadvantages of 
each method the reader is referred to an excellent dis- 
cussion by Dunham (1950). 

The refraction of energy in the lee of an island is 
difficult to determine by using a graphical approach. 

For this reason analytical solutions for the rays have 
been developed, where the depth variation is expressed 
as an analytic function (Arthur, 1951). Before proceed- 
ing with a discussion of this particular island case 
the relationship between wave height and energy will be 



-3 



reviewed . 

If it is assumed thac the total energy trans.uittec 
between two orthogonals is constant, convergence of the 
rays denotes a concentration of energy and a correspond- 
ing increase in wave heights. To a good approximation , 
wave height is proportional to the square root of the 
energy provided that the waves are not near the break- 
ing point (hiur.k and Traylor, 1947). The wave height;, 
however, will not usually be the same along a particular 
wave crest due to refraction effects over an irregular 
bottom. 

The mean wave energy per unic surface area for 
progressive waves is given by 

S = 1/8 pgH 2 

where p is the density of sea water, g the acceleration 
due to gravity and H the wave height. The mean wave 
power, EVb is constant between two wave rays and it follows 
that the relationship for the wave height is given by 

(1.3) H = K K K , 

r s 1 

where H, is the deeo water wave height, K is she 
refraction factor and X is the shoaling factor. 
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The refraction factor is defined by 

(1.4) K r = (b 1 /b) il 

where b is the distance between adjacent rays measured 
normal to the rays and the subscript 1 refers to deep 
water . 

The shoaling factor can be expressed by 

(1.5) K c = (V /V) 3 * 

where is the deep water group speed. In general 

V = nC 



where n is given by 



n 




4'.'V'h/L 
sinh 477’" h 



o- 



For the wave periods considered in this study n is 
essentially unity for all depths and hence V = C. Thus 

(1.6) K = (h./h)^ o 

s 1 

For an incident wave of unit amplitude it follows 
from (1.3) and (1.6) that the relationship for relative 
amplitude is given by 

A - K r (h 1 /h)" 4 . 



(1.7) 



15 



2. Fermat's Principle 

Fermat's principle of lease time is discussea in 
Joos (1934) o It can be statec thus: the propagation of 
light always takes place in such a way that the actual 
optical path (e.g. length of geometric path multiplied 
by index of refraction of the medium) is an extreme 
value compared with all other paths which do not follow 
the law of optics. In a wave refraction analogy, we allow 
the light rays to correspond to orthogonal s in a gravity 
wave refraction diagram. 

The application of Fermat's least time concept to 
refraction by concentric, circular, bottom contours is 
discussed by Arthur (1946). The travel time along a path 
from a point A to a point B of a refracted ray is given 
by 

r 

(2.1) 1=1 ds/C , 

J A 

where C = ( gh ) 2 along the path and ds is the arc length 
along the path. Consider depth h = h(r) only; then for 
I to be a minimum requires an r(-S-) such that <£l = C. 

For r>r^, the depth is assumed constant and che or- 
thcgonals are straight lines. For r <r<r. the r Q y s 
are curved. Using the coordinates as shown in Figure 



3, (2.1) takes the form 
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G. 



( 2 . 2 ) 



B 



I = 






dr / q-0-) ' 



2 _ 
r Qt; 



C(r) 



A 



A path which makes the integral of (2.2) a min- 
imum is found by the methods of uhe calculus of varia- 
tions (Joos, 1934), which in this case must satisfy 
the Euler-Lagrange condition 



i 5 ? 

d 1 dr/de' - r + G/C dC/dr = 0 , 
dG [ GC j GC 



where G ^ 
Let F 
only on r 




= G/C ( r ) . For C 
and r ' then 



and r ’ = dr/dG- . 
independent of -G- 



F depends 



(2.3) 



d_ 

dG- 



- r ' 



5>f 

Jr 1 



- 0 



r* / 0 , 



which implies that 



F _ r - 3X 
- Pr- 



= K, a constant 



If the constant of integration is evaluated at 
r = r^ the equation for the ray path takes the form 



(2.4) + dG 



dr/ 




2 (l/sin 2 G- 1 ) (h., /h(r ) ) 




3 



the sign being chosen according to the direction or the 
wave ray* 
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3. Kiv Path 



Consider the ray equation 



(3.1) 



de- 

dr 



r \ / I r -1 



where 



h 



2 2 
r, sin e- 1 



Now h = mr for r <rrcr, where m is the isxar.d 

o 1 ’ 



slope. Thus 



(3.2) 




o 



Substituting (3.2) into (3.1) and carrying out 
the integration, with the boundary condition -G- = -& 1 
at r = r 1 , gives 



(3.3) 



CSC 



■e-i + -G- 



- 2 _ 
r 1 sm -6- 



The wave ray which enters at r. , -3-^ will react, the 

island boundary, r , at -G- = -Q- q where 

f r „ ) 

, — — - sin ' 

(3.4) -e. = - e-, + 2c tn j r 1 \ , 

o 1 \ . ’ 



sin •©-. 



J 



_ o 



2 

provided that sin o 

2 

Consider sin -0% - r /r. , 

-7T--0-! 



■chen 



9 

For sin“ /r^ , the ray does not reach the island 

perimeter. In this case, the ray becomes tangent to a 
circle of radius r^ given by 

(3.5) r = r, sin^ -G-, , 

c l i 



i.e. r => r at -0- = - -€h. + 77* • A limiting ray exists ac 
co 1 ' 1 

r^ maximum or when -0- = 7T/'2 (see Figure 4 )„ Beyond this 

maximum critical radius, wave rays are unaffected by the 
island. 

For r = r , the inicial critical ancle can be decer- 

c o ’ 

mined by 

(3.6) sin -G-. 

lc 




For the particular island configuration considered here, 

= 20.44°. Hence rays will converge on the island 

oerimeter for -159.56°— ■©• X 159.56° . 

o 
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(4.1) 



4. Refraction Factor 

The refraction factor is given by the equation 

~ r dG cos Cl 
-2 o o ^o 



K 



r^G^cos G^ 



where r^cos G-^dG^ is the mathematical expression for the 

rav separation at r, and r dG co sCY represents the wave 
- ^ loo ^o 



seoaration at r . From (3.4) 
- o 



(4.2) 



dG 

c 

dG. 



2cos -G. 



1 



— - sin 2 G 
r l 1 



- 1 



From Figure 5. 



raG 



tan a = dr 
Substituting for dG/dr from (3.1), 



tan CC 



2 



n 



- 1 



By further trigonometric substitution for 

\ 35 



(4.3) 



cosCX Q = |l- 



. 2 . ^ \ ^ 

r i sin n 



/ 



Substituting (4.2) and (4.3) into (4.1) gives 



h 2 



r 



2 - sec 9^ *~\/ — — - sin 2 ^ 



This is valid for G-^ 



— G ~ <■ 

-A- G- 



20 



5. Phase Lag 



The water level displacement of a sinusoidal wave 
is given by 

(5.1) ^ = A cos (u/t - (j) ) 

where (J) is the phase lag at any point on the ray rela- 
tive to some fixed reference. The change in (J) along a 
ray can be expressed mathematically by 



(5.2) 



dd) _ UJ 
ds " C(r f 



where ds = 



o ■> 1 

2 >2 



(dr) -r (rdG-)' 

The phase change along the path from r = r. to 
r = r Q using (3.1) and (3.3) is given by 

2 L Ufa I — 5 



(5.3) 






— — • l/r.-r, sin -G-, - ~\ r -r, sin ■©-. 

1 • v j 



A correction f\ 0 must be added for the change be- 
tween the reference line <J) = 0 and r = r^ (see Figure 6 ) 
where 

tT r 1 (i - cos o 

Vg h ! 



(5.4) /\(j) 



Thus for the total phase lag as a function of -G-. 

j 
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( 5 . 5 ) 



(j) = 



2u/r. 



r 



V^7 L 



1 + cos -G--, ) 




. 2 
- sm 



1 

J 
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6. Relative Rower distribution 

The general relationship for the refraction factor 

for r < rcr. is given bv 
o 1 




.Substituting for r from (3.2) gives 

2 ^G.^ + -S-^j G^ + G- 

(6.1) K = sin -6 esc \ ■ — r — 2-tan -G- n ctn — - — g — , 
r 1 2 1 2 ■ ’ 

l U J 

valid for G-j <• G- <■ ( 2 7T - 3 ■©•-,) . for G = G. , K =1. For 

G = 'ff - G., , e.i. the critical angle for a minimum r of 

a ray, K“ 2 - 2 sin -G^. For G = 27T - 3 G^ , K~ 2 = 3. 

For direct rays characterized by a particular value 
_2 

of G. , is given by (6.1) for any -6- along the ray 

path. For those rays with |gJ <=■ G n , reflection occurs 

at r = r , = G q « The reflected ray passes through r = 

r. at G = 2 G- - G^ (see Figure 8). It emerges straight 

bevond r. and makes an angle 2 G from the x-axis. Thus 

two rays G 1 and G^ diverge at angular spread 2(G^ - G q ) 

and their intersection always lies within a circle of 

radius r. . Thus for r > > r. at |g| = 2 G- and Ig. <=• 

1 1 1 o 1 

1c 



I 
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P r, cos ©-, d©, 

K = — - - 

r 



Substituting from (4.1) 



2rd©- 

o 



( 6 . 2 ) 



K 2 = 
r 



2r 



r^cos -©i 



2cos ©-, 



VF 7 



- 1 



. 2 , 
sm ©-, 



For those rays with 77/2 ^(©-J >©-^ c no reflection 
occurs. The rays are refracted and emerge at r = r^, ©- 
= 2 Tf - 3 as a straight line and make an angle of 
2 Tf - 2 ©-^ , from the x-axis. Thus two neighboring rays 
©•^ and ©-^ have a divergent angle 2(-©\j^ - ©-^ ) and intersect 
within a circle of radius r^. For rs-^r^ at )-©-[ = 2Tf 



- 2 *1' *lc 



(6.3) 



I-©-} = 77/2. Then 



, r 1 cos 

r 2r 



For any given radial distance, r, from the center 
of the island, the relative power intensity can be 
expressed in the form 

(6.4) F (©•) = „2 

r l K r * 

The form of this "beam pattern" for reflected and refrac- 
ted waves in the far field is determined by (6.2) and 
(6.3), respectively. The total, relative power intensity 
in the far field is the sum of that contributed by the 
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reflected and refracted waves along raaial ~ines of -tr 
values (see Figures 9 and 10 in Appendix II). 
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CHAPTER III 
DISCUSSION OF RESULTS 

A calculation summary for the circular island 
aopears in Appendix I. This is presented in two parts: 
(A) information pertinent to the rays and refraction 
on the island lee shore; (3) information pertinent to 
the far field power distribution. The graphical repre- 
sentation of phase lag, refraction factor, relative 
amplitude, and relative oower intensity for the island 
shore are given in Figures 9 to 12 (Appendix II), A 
simple wave refraction pattern consisting of direct 
rays into the island shore is produced by the constant 
slooa case. For those incident waves entering with 
azimuths within the range from -20,4 to 20,4 degrees, 
the waves converge on the island covering an angular 
range from 0 to - 159,56 degrees, A wave shadow zone 
of about 40 degrees range is formed in the lee of the 
island and is an area unaffected by refraction effects. 
The shadow zone is developed between the shore and the 
outgoing critical ray. 

"or those rays with j -©^ 1 <• -©-^ , reflection occurs 

at r = r .-©- = -©• and the ray emerges from r n at •©- = 
o o 1 

2 -e- - -9-, , making an an r ’e of 2 -9- with the x-axis. 

o i. o 

For those rays with 7T/2 >9^, no reflection 

occurs, but the rays are refracted and finally emerge 
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at r = r\ , -&■ = 2 77" ~ 3 -G-^ and with an angle of 2 7 7 
-2-G-. from the x-axis. At cha ooint of emergence K 
= l/~\/~3, however, at -©^ = - ']~{'/2 v a discontinuity exists. 
Incident rays that are just outside r n at - 77"/^ are 
unaffected by the island, while those rays just inside 
r n undergo refraction. 

_L 

For those emergent rays at large distance from one 
island (r> >r^), the center of the island can be con- 
sidered as the origin of these rays and the relative 

2 

power intensity in the far field (rK /r n ) is simply a 
function of •©-. 

The relative power intensity is found to increase 
gradually from [■©•} = 180 to an extreme at 2 -S-. (40.8 

degrees) at which there is an abrupt decrease. This 
discontinuity in the far field beam pattern of the 
power is caused by the shadow phenomenon inherent in 
this refraction analysis. The contribution of the re- 
fracted waves in the far field exceeds that of the re- 
flected rays by more than a factor of 4 (see Figures 13 
and 14). 

On the other hand, when considering the wave 
amplitude at the island shore, derived by use of Green’s 
Law, the shoaling factor is significant while refraction 
effects are limited by the small -G-^ value. Furthermore, 
for a constant slope of 0.1, any variation in the 
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island contour parameters will not appreciably c._ ter 
the critical angle. The refraction factor varied a srr.ux_ 
amount due to the slope and parameters of the island 
selected (see Figure 10). 

The variation in relative amplitude, in this case, 
is governed by the refraction factor and drops off 
abruptly to zero at the point of ray tangency to the 
island. The energy for the near critical rays coes not 
dissipate at a single point on the island shore as 
suggested by the refraction analysis. It is in this 
region of energy convergence that the assumption of 
constant power between wave rays is no longer valid 
and a down gradient flow of energy occurs across ort’nog- 
onals. The "shadow zone" of the island is no longer a 
region of constant energy flux, but would contain energy 
input lost by diffraction effects near the point of 
tangency of the critical ray. Thus, by conducting a re- 
fraction analysis alone, it is difficult to adequate- 
ly describe the energy distribution on the wave lee 
side of the island. This also applies to the effective 
shadow zone in the far field pattern. 

On the far side of the island, where a diverging 
cross-over pattern of emerging rays is formed, the 
relative amplitude can be derived by combining that 
of two outgoing wave rays at tneir proper phase relation- 
ship. 



with 



The phase lag on the island shore increases 
- 3 - q (see Figure 11). The range of phase lag at the 
island shore varies from about 1 radian for v/aves of 
12 minute period to about 6.5 radians for waves of 2 
minute period. Thus the phase lag of waves in the lee 
of the island varies considerably with period, while 
one amplitude is unaffected by the period (within the 
tsunami range). 
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